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Motivated by recent development in quantum fidelity and fidelity susceptibility, we study relations 
among Lie algebra, fidelity susceptibility and quantum phase transition for two-state system and 
the Lipkin-Meshkov-Glick model. We get the fidelity susceptibility for SU(2) and SU(1, 1) algebraic 
structure models. From this relation, the validity of the fidelity susceptibility to signal for the quantum 
phase transition is also verified in these two systems. At the same time, we obtain the geometric 
phase in these two systems in the process of calculating the fidelity susceptibility. In addition, the 
new method of calculating fidelity susceptibility has been applied to explore the two-dimensional XXZ 
model and the Bose-Einstein condensate (BEC). 



PACS: 03.65.Fd, 03.67.-a, 03.65.Vf 

1. Introduction 

Fidelity, one of the most intriguing feature of quantum information science,' 1 ^ has been 

widely studied in recent years. [ 2 ~ 14 ^ The quantum phase transition(QPT), driven purely by 

quantum fluctuations and occurring at zero temperature, is believed to be an important concept 

in condensed physics. ^ In a quantum many-body system, the QPT is driven by purely the 

quantum fluctuation in ground states. The wave function of the ground state can have a 

abrupt change as the system varies across the phase transition point. Therefore, an approach 

to quantum phase transitions based on the quantum-information concept of fidelity has been 

put forward. ^ However, except for a few specific modelsj 2 ' 16 ' the calculation of the ground-state 
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fidelity is tedious. Recently, a neater and simpler formalism of fidelity to critical phenomena 
was introduced, ^ for the so-called fidelity susceptibility, to signal the whole QPT. The main 
advantage of this approach lies in the fact that the fidelity is a purely Hilbert-space geometrical 
quantity and no a priori knowledge of the structure of the considered system is required for 
its use. As shown in Ref. [17], the fidelity susceptibility is intrinsically related to the dynamic 
structure factor of the driving Hamiltonian that is evaluated though the scheme based on some 
numerical techniques including exact diagonalization and density matrix renormalization group. 
On the other hand Zhang et alJ 18 ^ employed Lie algebra to evaluate the fidelity susceptibility, 
and show high efficiency of it. 

In this paper, we investigate, in a general framework, how the differential form of the fidelity 
susceptibility can be established in terms of the general Lie algebra. As will be seen below, 
under certain conditions , the differential form is an effective tool in detecting the critical 
points of the QPT. To demonstrate this, analytic formulas for the fidelity and the fidelity 
susceptibility are derived for the SU(2) and SU(1, 1). By using these formulas, the fidelity and 
the fidelity susceptibility can be easily calculated for a large class of many-body systems, as 
long as the Hamiltonian of the system can be rewritten as the form [see Eq. (4)]. Employing 
our general formulas to the two-state system, we can show the fidelity susceptibility of the 
system in terms of SU(2). At the same time, the geometric phase, which is also an effective 
indicator in detecting the QPTj 19 ~ 22 l can be obtained in the process of calculating the fidelity 
susceptibility. On the other hand, according to the general expressions of the differential form 
of the fidelity susceptibility, one can also expect that the same results for the Lipkin-Meshkov- 
Glick(LMG) model in terms of SU(1, 1). Furthermore, we extend this differential form fidelity 
susceptibility to other physics models. 



where Hj is the driving Hamiltonian and A denotes its strength. The fidelity is the absolute 
value of the overlap between two ground states |\Po(A)) and (^(A + <5A)) 



2. Formulism 



The general Hamiltonian of quantum many-body systems reads 



if (A) = H + XH r , 



(1) 




(2) 



with 5X a small deviation. Extracting the leading term of the fidelity, the fidelity susceptibility 
can be obtained' 17 ' 

with eigenvalues E n (X) and corresponding normalized eigenvectors |\& n (A)). The eigenstates 
define a set of orthogonal complete bases in the Hilbert space. 
Given a physical system, the Hamiltonian can be written as 

H = Y,£iHi + E( A «£« + KE- a ), (4) 

i a 

where and X a are coupling parameters and {Hi, E a , E_ a } are the Cartan-Weyl basis of a 
semi-simple Lie algebra. Using the unitary operator U(£ a (X)) = exp[J2(^ a E a — £*E- a )], the 

a 

Hamiltonian can be turned into the diagonal form 

&t(Z a )HU(S a ) = J ErjiH i . (5) 

i 

Therefore, the eigenstates of Hamiltonian (j3J) are \^f) = U(£ a )\ref), where \ref) is the direct 
product of normalized eigenstates of Hi. Then the absolute value of the overlap between 
|*(A)) = U(C{X))\ref) and |*(A + 5X)) = £>(£ a (A + SX))\ref) is 

F=\(ref\U\U^)U(U^ + SX)))\ref}\. (6) 

According to the general relation between fidelity and fidelity susceptibility, F = 1 — |(<5A) 2 x(A) + 
• - -, the expression of the fidelity susceptibility is given by! 18 ! 

X (X) = -(ref\(U^d x U) 2 \ref) - \(ref\Wd x U\ref)\ 2 . (7) 



3. Two-state systems 

The two-state system is the simplest quantum system which can be calculated exactly. 
Furthermore, the two-state systems possess several typical quantum properties. So we shall 
study it firstly. The Hamiltonian of a two-state system can be written as 

H = — B • a, (8) 

where B is an external magnetic field and a = (cr x , a y , a z ) are Pauli matrices. In a z basis 
| f), | l), Pauli matrices take the form 
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We can rewrite H in the SU(2) form 

H = -2Bcos6J z - Bsin6e^J + - BsinOe-^J-, (10) 

by means of the generators of the algebra SU(2) 

111 

J z = -&z, J+ = -cr+, J- = 2 a ~- I 11 ) 

These satisfy the usual commutation relations 

[J+, J_] =2J Z} [J Z ,J ± ] =±J±. (12) 

Resorting to the Eq. (J5J, we introduce the unitary operator 17 (6 1 , 0) = exp(— |e~^ J + + |e l9i J_). 
The Hamiltonian (JSJ) then can be diagonalized 

1^(0, (f>)HU(0, 0) = -5J Z . (13) 

Here, and can be regarded as adiabatic parameters. For simplicity and without loss of 
generality, we fixed first and have 

U\6,(j))d <f> U{e,(f)) = l -sm6{e-^J + + e^J_) + 2*sin 2 9 -J z (14) 

and 

| {ref\U\e, 0)^*7(0, 0) |re/)| 2 = 4 sin 2 ? (15) 

where \ref) is the cr 2 basis | f), | J,) i n the two-state system. Finally we get the fidelity 
susceptibility of spin-| subjected to an external magnetic field 

XF = \sm 2 6. (16) 

Since {ref\U\9 ) (f))d ( f ) U(9 ) (f))\ref) is just the Berry adiabatic connection, which contribute 
a Pancharatnam-Berry phase' 23 ' 24 ' to the spin as the magnetic field rotates adiabatically around 
cone direction, we get the Berry phase 

7(M) = ~i fj(ref\U\e^)d^U{e^)\ref) 

= ±tt(1-cos0). (17) 

The results of Eqs. ( TTBl) and ( JT7|) are in agreement with that in Ref. [25], which verify the 
reliability of our method. 



4. The Lipkin-Meshkov-Glick model 

The LMG model' 26-28 ] was originally introduced in nuclear physics. It provides a simple 
description of the tunneling of bosons between two degenerate levels and can thus be used to 
describe many physical systems, such as two-mode Bose-Einstein condensates' 29 ' and Josephson 
junctions J 30 ' Recently the entanglement in this model has attracted much interest because of 
available numerical calculations and plentiful phase diagrams. ' 31 ' In the thermodynamic limit, 
its phase diagram can be simply established by a semiclassical approach J 32 ' For finite large 
number N of particles, it was studied by the 1/N expansion in the Holstein-Primakoff single 
boson representation' 33 ' and by the continuous unitary transformation.' 34 ' 

The Hamiltonian of the LMG model can be written as 

i<j i 

= -^(S2 + 7Sj)-2/»S, + ^(l + 7 ) 
= -^(l + 7)(S 2 -52-y)-2/i5, 

-2^(1-7)051 + ^), ( 18 ) 

where the <j k {k = x, y, z) are the Pauli matrices, S K = J2i c«/2, and S± = S x ±iS y . A and h are 
the spin-spin interaction strength and the effective external field, respectively. N is the total 
spin numbers and 1/N ensures that the free energy per spin is finite in the thermodynamical 
limit. It is understood that H preserve the magnitude of the total spin and the parity P = Hi & l z , 
i.e., 

[H,S 2 ] = 0,[H,P]=0, (19) 

for all values of the anisotropy parameter 7. Specially, in the isotropic case 7 = 1, one has 
[H, S z ] = 0, so that H is diagonal in the eigenbasis of S 2 and S z . We adopt the 1/N expan- 
sion method corresponding to the large iV limit. We first use the Holstein-Primakoff boson 
representation of the spin operator' 33 ' in the S = N/2 subspace given by 

c c t N t 
b z = b — a [ a — — — a'a, 

S + = (2S-Ja) 1 *a = N^l-^) 1 ia = SU (20) 

where the standard bosonic creation and annihilation operator satisfy [a,a^] = 1. This repre- 
sentation is well adapted to the computation of the low-energy physics with (a)a)/N 1. 



The next step consists in inserting these expression in Eq. ([18]) . and to expand the argument 
of the square roots. Keeping terms of order (1/iV) -1 , (l/iV)^ 1 / 2 , and (1/N)° in the Hamiltonian 
yields (h > 1) 

H = —hN + (2h-l- 7 )a f a - [(1 - 7) /2](a t2 + a 2 ). (21) 

In order to analyze the fidelity susceptibility of the above equation, we transform Eq. ( 12TT) 
into the form of SU(1, 1) following the ideas in Ref. [35]. We first introduce the generators of 
SU(1,1), 

K+ = = l -a 2 ,K z = ^(2aW 1), (22) 

which satisfy the communication relations of SU(1, 1) algebra. Submitting these expressions of 
the SU(1, 1) generators into Eq. (l2Tj) . one can get 

H = -hN - Uh - 7 - 1) + 2(2/i - 7 - l)K z + (7 - l)K + + (7 - 1)K_, (23) 

which consists with Eq. (jl]). We do not show the diagonalization of SU(1, 1) algebraic structure 
model explicitly here, but interested readers are recommended to refer to Ref. [18]. We simply 
quote the main result that is connection with the fidelity susceptibility of the SU(1, 1) algebra 
model. After all the procedures, the fidelity susceptibility of SU(1, 1) algebraic structure model 
becomes 

X^ (1) i) = g[(f) 2 + sinh 2 ^(|)]. (24) 
Assuming tanh^ = (1 — j)/(2h — 1 — 7), the fidelity susceptibility is represented by 

XF ~ 32(l-/i) 2 (/ i - 7 ) 2 ' 1 ' 

which is in agreement with Ref. [36]. 



The derivation above is only valid for h > 1, for < h < 1 the calculation is actually similar 
to the above case of h > 1. When < h < 1, tanh# = 2 ^ fa l2 7 > the fidelity susceptibility is then 
obtained accordingly 

= ^rh^F (26) 
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Fig.l. Fidelity susceptibility as a function of h and 7. The divergent character of xf is 
clearly displayed as h — > 1. 
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Fig. 2. Fidelity susceptibility as a function of h for 7 = 0.5 in large N limit. 

Thus we obtained fidelity susceptibility of the anisotropic LMG model in the large N limit. 
As shown in Figs. 1 and 2, it is obvious that xf is divergent at h — 1 where the LMG model 
has been proved to experience a second-order phase transition independent of the anisotropy 
7, which is well described by a mean-field approach.^ Through the calculation of fidelity 
susceptibility, we have found the geometric phase. For illustration, we consider the system which 
has a rotation g((p) around the new z direction. The Hamiltonian becomes H((f)) = g((p)Hg\(j)) 
with g(<j>) = e^ s *. Then Eq. (J2D can be written as 

H = -hN + (2h-l- 7 )a f a - [(1 - 7)/2](a t2 e" 2 ^ + a 2 e 2 ^). (27) 

The geometric phase of the ground state, accumulated by varying the angle <fi from to it, is 
described by /3 = — i J$ {g\-§z\g) Finally, we get the geometric phase 

P =~i £(9\§i\9) = -1 £(1,0|C/(M)WM)|1,0> 

= -i [ n -2ismh 2 - = vr(l - coshfl). (28) 
Jo 2 
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The derivation above is also only valid for h > I. When < h < 1, tanh# = 2 fc ^2 7 ■ Therefore 
we can get geometric phase of the LMG model in the thermodynamic limit. One can find f3 g is 
also divergent at h — 1 which is equal to the fidelity susceptibility in indicating the quantum 
phase transition in Figs. 3 and 4. 




Fig.3. The geometric phase as a function of h and 7. It is obvious that /3 g is divergent at 
h = 1 independent of 7. 




Fig.4. The geometric phase as a function of h. The parameter 7 = 0.5. 



5. Other models 



In this section, we generalize the new method of calculating fidelity susceptibility to other 
quantum many-body models. As there is nobody to obtain the analytical results of the fidelity 



susceptibility of the two-dimensional XXZ model and the Bose-Einstein condensate precisely, 
we expect the under results can arouse more wonderful ideas or results. 



5.1. Two-dimensional XXZ model 

The Hamiltonian of the XXZ antiferromagnetic model reads 

H xxz = £(Sf Sf + SfS? + riS?S?), (29) 

(ij) 

where Sf(a = x,y,z) are the spin-1/2 operators at site i and rj = J Z /J X {J X = J y ) is a dimen- 
sionless parameter characterizing the anisotropy of the model. The sum runs over all the nearest 
neighbors on a square lattice. We begin with the two-sublattice model and Holstein-Primakoff 
transformation, 1 33 1 

s: = V2s a t(i - ^) 1/2 ,^- = (s:)\ 

S z a = S + a f a, S z h =S- tfb, (30) 

where a\ a{b\b) are boson creation and annihilation operators on sublattice A (sublattice B), 
respectively. The particle numbers n a = a^a, = b^b cannot excel 2S*. Transforming the 
operators into momentum space, we rewrite the Hamiltonian as 

Hxxz = -2 V zNS 2 + 2zSJ2Hk, (31) 

k 

where z is the coordination number of the lattice and Hk is of the form 

H k = v(aia k + &t&k) + 7kK&k + 4 & k) • (32) 
Here 7^ = z~ 1 J2 etkn with R is a vector connecting an atom with its nearest neighbor. We 

R 

can rewrite if k in the SU(1, 1) form, 

H k = 2r ] A^ + lk (Al + A k ), (33) 
by means of the generators of the algebra SU(1, 1) 

-4+ = 4&t, A* = a k b k , 

A* = \{K + 4 + 1). (34) 



Taking a similar transformation of H^, we get the unitary operator U = exp(£k^4+ — 6c^._) with 
tan^k = — -jjj 1 - With the spin-wave theory framework, we obtain the unitary operator and then 
calculate the fidelity susceptibility of the model in two dimensions. Taking tanh# = e~ l * = 
2k = 1 into the Eq. (|24|) . one can obtain 

f {ydxlk - Ikdxrj) 2 dk 
X J 32(^- 7 ^ (2vr)3- ^ 

5.2. Fidelity susceptibility for the Bose-Einstein condensate 

The standard description of the Bose-Einstein condensate is via an order parameter field 
ty(x). The Hamiltonian takes the standard form 

r h 2 

= j d 3 x[— \W(x)\ 2 + U(x)\V(x)\ 2 } 

+ \S d * x I d3 y^*(y)^*( x ) v ( x ^y)^(y)^( x )- ( 36 ) 

As in Ref. [35], the Eq. fl36|) can be rewritten in the form of SU(1, 1) algebra 

n = 2[a A 3 + ^(u A + + u*A°_)] 
+ EMs + \^ k A% + < A - )] " E*> (37) 
where a = e + | E (V + V r fc )((ra fc ) + (n_ fc )), m = V (ag)+ E Vfc(a fe a_ fc ), a fc = e k + (n )(V + V k ), 

k^O k^O 

u k = Vfc(ag), E* = UV \(a 2 }\ 2 + a ] + \ E [(cr* - e*)( n fc + + 0fc] + | E (u k (at a -k) + 

k^O k^O 

u* k (a k a_ k )). Introducing the generators of the algebra SU(1, 1) 

Al = l(n + \),Al = ^-,A°_ = | (38) 

and 

A t = \i n k + n_ fc + 1), 4* = a£aj; fc , = a k a_ k . (39) 

We calculate the fidelity susceptibility of Ho = 2[aoA 3 + \ {uqA° + + UqA°_ )]. From this equation, 
we can get 

[/ = exp(£ A° + - f M° ), £ = r exp(z0). (40) 

Taking the same procedures as in the two-dimensional XXZ model, we can get the fidelity 
susceptibility of the %q 
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Similarly the fidelity susceptibility of the T-Lk can be obtained 
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6. Conclusions 

In conclusion, we have established the differential form of the fidelity susceptibility in terms 
of the general Lie algebras. Meanwhile we investigate the geometric phase which can also 
show the phase transition point. Therefore, we construct the relation between the fidelity 
susceptibility and geometric phase. We also apply the differential form of fidelity susceptibility 
to other physics models. In particular, we focus on the SU(2) and SU(1, 1) algebras. The form 
of the fidelity susceptibility of SU{2) and SU(1, 1) algebra not only enables us to evaluate the 
fidelity susceptibility easily, but also builds a straightforward connection between quantum- 
information theory and the Lie algebra in quantum many-body physics. 
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